ASYMPTOTIC FITNESS DISTRIBUTION IN THE BAK-SNEPPEN MODEL OF 
BIOLOGICAL EVOLUTION WITH FOUR SPECIES 



ECKHARD SCHLEMM 

Abstract. We suggest a new method to compute the asymptotic fitness distribution in the Bak-Sneppen model 
of biological evolution. As applications we derive the full asymptotic distribution in the four-species model, 
and give an explicit linear recurrence relation for a set of coefficients determining the asymptotic distribution 
in the five-species model. 



1 . Introduction and main result 

The Bak-Sneppen (BS) model of biological evolution has been introduced in [2] to study self-organised 
critical behaviour in populations with natural selection and spatial interactions. In the model, a population 
is spread out over a circle so that each species has exactly two neighbours, and each site is assigned a 
numerical parameter between zero and one that describes the fitness of the species at that particular site. 
The system evolves according to the following discrete-time dynamics: at each time step the species with 
the smallest fitness parameters becomes extinct and is replaced by a new species whose fitness parameter 
becomes independent uniformly distributed on the interval [0,1]. To take into account interactions between 
species, e. g. competition for resources or predator-prey relations, the fitness parameters of the two neigh- 
bours of the least fit species are also reset to random values. This represents the idea that the fitness is not 
just a fixed property of a species but also depends on its environment. Even though this is a very rudimen- 
tary model for evolution it does possess some of the characteristic features thought to be associated with 
evolving biological systems, including long-range dependence and evolutionary activity on all time scales, 
see e. g. [10] for a discussion in the context of self-organised criticality. 

Despite its apparently easy definition, the BS model has withstood most attempts at mathematical anal- 
ysis in the past. Partial results have been obtained, however, in the context of rank-driven processes and 
mean-field approximations [3, 7, 8]. Early on it was conjectured based on simulations that the fitness dis- 
tribution at a fixed site converges, in the limit of large populations, to a uniform distribution on the interval 
[f c ,ll where / c * 0.667. 

There are only few rigorous results about the original BS model; in [12] it is proved that the expected 
fitness at a fixed site in the stationary regime is bounded away from unity uniformly in the size of the 
population; [13] gives a conditional characterization of the limiting distribution in terms of a set of critical 
thresholds [see also 5, 6]. 

In this paper we suggest a novel approach to determining the stationary fitness distribution in the BS 
model with a finite number of species. The method consists of translating the evolutionary dynamics 
into a set of linear recurrence equations for the coefficients of the densities of the finite-time distributions 
and analysing the asymptotic properties of these recursive equations. We apply the method to derive the 
asymptotic fitness distribution in the BS model with four species, which is given by (Theorem 1) 

It has been suggested that a similar, and potentially easier, analysis can be applied to the anisotropic 
Bak-Sneppen (aBS) model; in this variation of the model the evolutionary dynamics are modified so that 
the fitness parameter of only one neighbour, say the right one, is updated at each time step together with 
the fitness parameter of the least fit species. While many qualitative properties are believed to be shared by 
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the original and the anisotropic BS, analytical considerations and numerical calculations indicate that they 
fall in different universality classes, and have thus different critical exponents [9, 11]. Our method can be 
applied to the aBS. In the first non-trivial case of three species the asymptotic fitness distribution follows 
from computations which are very similar to those leading to Theorem 1, and is given by 

The complexity of the calculations for larger populations seems to be increasing at a similar rate as in BS 
and studying aBS does therefore not, for the purpose of understanding the method we propose in the current 
paper, offer major advantages over the original model. For this reason we concentrate in the following on 
the isotropic case. 

Outline of the paper. The paper is organised as follows. In Section 2 we present one of several possible, 
precise definitions of the BS model in terms of a Markov chain, and use this representation to derive 
qualitative properties for arbitrary population sizes. We then turn to the special case of the four-species 
model in Section 3, where we derive the full asymptotic fitness distribution. Finally, Section 4 is devoted 
to making further progress towards a similarly explicit solution for larger populations. In particular, we 
obtain explicit linear recursions for the coefficients that describe the finite-time distributions in the five- 
species model. 

Notation. Throughout the paper we denote by N the set of positive integers and write n e N for the number 
of species in the BS model. Arithmetic operations in indices pertaining to the number of species are always 
understood to be performed modulo n, so as to give a result between 1 and n. The symmetric group on n 
letters is denoted by S n . For a vector x e [0, 1]" and an integer v e {1, ...,«} we denote by (•)[>] and (-)]v[ the 
inner and outer v-centred projections, that is 

0[v] : [0,1]" -> [0,1] 3 , x [y ] = {x v -\, x v ,x v +&, 
(•)]v[ : [0,1]" -> [0,1]"" 3 , x M = (x Y+2 ,...,x n ,xu...,x v -2). 

The /th components of the vectors X[ v ] and Xj V [ are denoted by x^y and X] V y, respectively. Given another 
vector £ e [0, l] 3 we extend the notation of outer projections by setting 

(•)]v[ f : [0,1]" -> [0,1]", x ]v[f = (tji,&&,x v+ 2,---,x n ,xi,...,x v -2). 

We use the multi-index notation x', x e [0, 1]", i £ Nq, as a shorthand for the monomial x'^ ■ ■ ■ x' t " . We further 
let P denote probability and write SS(X) for the Borel sets of a topological space X. Finally, we denote by 
the indicator of an expression S, which is defined to be one if £ is true, and zero else; for a set S , we 
also write l$(x) instead of 1{ x <eS)- 

2. Definition of the model and qualitative analysis for general population sizes 

We interpret the Bak-Sneppen model with n species as a Markov chain f n = (/„ ^)a>o on the hypercube 
[0, 1]", see, e. g. [14] for a general introduction to this topic. To avoid the trivial case in which all fitness 
parameters are updated at each time step, we assume that n > 4. For each time k e N, the values of 
fnk = (fn,k,v) v =i n denote the fitnesses of the different species in the population. In order to formalise 
the heuristic description of the evolutionary processes described in the introduction we first observe that 
the fitness landscape /„ k+l at time k + 1 can be determined easily if one knows the fitness values at time k 
and, in particular, which species was the least fit at that time. In fact, if at time k the vth species is the least 
fit, then the fitness values of the species at the (v- l)th, vth, and (v+ l)th site are independent uniformly 
distributed at time k+l, while all other fitness values remain unchanged. For this argument to be valid, 
it is necessary that /„ k has a unique minimal value. In fact the BS model, as we have described it, is not 
well-defined without a rule for breaking ties. One natural possibility is to randomly select one of the sites 
with minimal fitness value. We will, however, not have to concern ourselves with this complication because 
we will assume that the initial fitness distribution at time k — is absolutely continuous, in which case ties 
occur with probability zero. Further formalising the argument, we can partition the event [f n i + \ e A), 
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A = n v [J«v".4] e #([°. !]")' as U"=i {/«,fc+i e A) n{min/„ >jt = f nXv ] and obtain, 
P(/„,*+i e a) =X P (^+i 6 A and nfofn* = 

v=l 

= 2 P (Wm 6 Ajvl and = /"Av)P((t/i,I^2, f/3) e A M ), (2.1) 



v=l 



where Q] v [ and (•)[>] denote the inner and outer projections, respectively. The random variables U\,U2,Ui 
are independent and uniformly distributed on [0,1], and thus the last probability equals Leb 3 (A[ v j). Put 
differently, the transition kernels P„ x (-) of the Markov chain /„, which are characterized by the equations 

v(fnMi &A \fn,k = *)= f<*(d n £), *e[0,l]", Ae#([0,l]»), 
Ja 

are given by 

P^(d"^)= [] ^(d^)d 3 (^_i,^ v ,f v+1 ), v = argrninx; (2.2) 

/i£{v-l,y,v+l} 

here, (£,*) e [0, 1]" x [0, 1]", and 6 X (-) denotes the Dirac measure located at x. This quite explicit expression 
is enough to prove that the Markov chains f n are ergodic in a strong sense, which is crucial for our approach 
to determining their stationary distributions. For the definition of uniform convergence we refer the reader 
to [14, Definition (16.6)]. 

Proposition 1. For every positive integer n, the Markov chain f n is uniformly ergodic. 

Proof. It is straightforward, albeit tedious, to show that the multi-step transition kernels ¥% x (-), which are 
defined recursively by 



JfeA Jve[0,11" 



dominate, for each m 3= n/3 - 1, a constant times Lebesgue measure restricted to [0, 1]". The claim is then a 
consequence of [1] [Theorem B, ii)]. Since a similar result has been obtained in [4], we omit the technical 
details. □ 

Proposition 1 implies in particular that there exists a unique invariant distribution /j„ satisfying 
H„(A)= f f P^(d n £K(d"*), forallAe^([0,l]"); 

JfeA Jx€[0,l]" 

moreover, the finite-time distributions, irrespective of the initial distribution at time k — 0, converge to fi n in 
the total variation norm. One way to obtain results about the stationary distribution p n is thus to fix a conve- 
nient initial distribution for /„ , to determine the finite-time distributions, and to analyse their asymptotic 
behaviour. The most convenient choice of initial distribution turns out to be the uniform distribution f U^ l y 
We now proceed to give a qualitative description of the finite-time distributions. 

Proposition 2. For every n > 4, the finite-time distributions of the Markov chain f n are absolutely contin- 
uous with piecewise polynomial Lebesgue densities. More precisely, for every k € N and o~ e S n -3, there 
exist finite index sets I n ,k,cr c Nq~ 3 and polynomials q n x<r E Q[*i , • • • , -*n-3] of the form 

q n ,kAx)= Yj a "Mx\ *e[0,l]"- 3 , (2.3) 

such that 

g n , k (x)d"x € d"*|/„ t0 ~ <U» tl ) 

n 

\d n x. (2.4) 

v=l o-e5„_3 
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Proof. The finite-time distributions are absolutely continuous because the initial distribution at time k — 
is, and the transition kernel 2.2 preserves absolute continuity. A reformulation of Eq. (2.1) shows that the 
densities g, hk satisfy the recursion 



gn,k+l( x ) ~ 2_i L , %2<™n(fi^3 
v=l ^[0.1] 



(2.5) 



The claim thus follows by induction on fe, using the fact that integrals of polynomial expressions over 
bounded convex polytopes are again polynomial expressions. More precisely, we assume that, for some 
k > 0, the function x h-> g n ,k(x) is of the form asserted in Eq. (2.4). It then follows that 



v,/i,cr,( "[0,1] 



• 2 <min(f 1 ,f 3 ^:] v[ )} J 



1/4 



(*]v[f) 



where, for j> e [0, 1]"~ 3 and cr e S„-3, the expression^ ~ cr is a shorthand for <y<r(i) < . . . <3'o-(«-3). Since, 
for each /j = 1, . . . ,n, the integrand is equal to /?(*] y [) for some polynomial function e Q[xi,., . ,JC n -3L the 
claim follows. □ 

For the special cases n = 4,5 we will see later how Eq. (2.5) can be translated into explicit recursions 
for the coefficients a, h k,cr,i- Before we do that, however, we compute, for general n, the density g n \ of the 
fitness distribution after the first evolution step. 

Lemma 1. For each n > 4, the density g„ t \ is given by 

g„,i(x) = l[o,i]"(x)^p(minx] V [), p(x) = x(^x 2 - x+ lj. 



(2.6) 



Proof. The claim follows from Eq. (2.5) and the observation that an initial distribution of 'W^y corre- 
sponds to g n fi(x) = 1 . Alternatively, one can use the fact that 

P(t/ 2 < min(C7i, U 3 ,y)) = |p(min(t7i, U 2 , U 3 ) <y) = i [l - (1 -y) 3 ] . □ 

3. Derivation of the asymptotic fitness distribution for four species 

In this section we expand on the results from the last section to compute the stationary distribution in 
the Bak-Sneppen model with four species, n = 4. This is much easier than the general situation and we are 
able to obtain a complete picture of the joint asymptotic behaviour of the four fitness-parameters. Partial 
generalisations to larger populations are discussed in Section 4. For notational convenience we suppress 
the index n — 4 in this section. Proposition 2 then implies that there are rational numbers a k ,i and positive 
integers dk such that the finite -time distributions of / = f 4 are given by 

[ 4 ] <k 

p(/ i ed 4 x|/ ~'Z/® 4 1] )= J]q k (x v ) l [0jl]4 (*)d 4 *, q k (x) = Y J a k ,,x i , 



Lv=l 



i=0 



so that gk(x) = g4,k(x) = TiiTiv=\ a k.' x 'v The Chapman-Kolmogorov-type equation 2.5 implies that, for 
k>0, 

4 4 r 

Y^Yj^^^YjYj"^ I , 1 fe<™n(f,,ft,^)l[4+^2 + ^i + ^] d3 ^ 
y=l i v=l i J[0 ' l] 

Equating each of the four v-summands separately and evaluating the integrals on the right side, one sees 
that, for v= L....4, 



s 



atk+uXy 



= V a k j x\ I ljf 2 <min(f,,ft,^)}d 3 ^+ I l {f2 < m in(^ 1 , ft ,A V )}^2 d3 ^ + 2 \ X \h 
i I J[0,1] 3 J[0,1] 3 v7[0,l] 3 

t-r 1 [\3 / \i + 3 i + 2 i + l) i + l\; + 3 ' + 2/J t-r 1 i + 1 
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Comparing coefficients of equal powers of x v on both sides of this equation, one arrives at the following 
recursion for the coefficients a^j. 

Proposition 3. For each k > 1, the constant term akfi vanishes. The remaining coefficients a^j satisfy the 
recursion 

-■3ir+l a Kj 



ask+l,i 



£ ^j=l j+l> l ~ L > 

-2^,1 + 3^,2, i = 3, 

'+1 ™ £+1 ™ i 1 i+l ; \ A 

— ak,i- 1 - jzj ^,1-2 + 3 j^a*,i-3 . ' > 4 , 



/n particular, a^j is equal to zero ifi^3k+l. 

Proof. Only the last statement requires further justification. By the last line of the recursion 3.2, the 
coefficient a^+ij can be written, for any £= l,...,k+l, as a linear combination of the coefficients cek+1-tj-m.y 
m = €,..., 3(. Choosing I — k + 1 and observing that ao, i is zero for i > 1, it follows that 0%+ 1 vanishes for 
i>3(fc+l) + l. □ 

The next step consists in finding as explicit a solution to this recursive system as possible. Fortunately, 
since we are mainly interested in the asymptotic behaviour of the Bak-Sneppen model, it suffices to con- 
centrate on large values of k. 

Lemma 2. For every i e N, there exists a rational number such that a^t = Pi for every k > In 
particular, lim^co a^i — Pi- 
Proof The claim is true for i = 1; this follows from the fact that gk : x i-» Yjv1k(x Y ) is a probability density 
function and thus integrates to one; by symmetry, qi< then integrates to 1/4. Observing that + 
equals J qk(x)dx, it follows from Eq. (3.2) that, for all k > 1, 

3/t+l , , 

~ V ak -i n 1 1 

j=i J 

Similarly, = 1/2, k^ 3, and a^i = -1/3, 4, and so the claim is true for i = 2,3 as well. For general 
i > 4, the claim follows by induction. □ 

Proposition 4. Let Pi = a^ij fee fne numbers from Lemma 2. The functions qt converge uniformly on [0, 1], 
i. e. in f/ze supremum norm, to q :— x i-> YtiPi x '> which is given by 

3 jc(3-x(3-x)) 
2(3-x(3-x(3-x))) 2 

Proof. It follows from Proposition 1 that the functions ^ converge uniformly to some limit, and because 
a/cj converges, for each i, to by Lemma 2, this limit must beiw YiiPi x> ■ Using Proposition 3 one sees 
that the sequence c, = (i + 1) arj+i.j satisfies the homogeneous four-term recursion 

Ci = Ci- 1 - c { -2 + - Ci-3 , i > 4. 
By the theory of homogeneous linear recurrence equations, c; is given by 

Cj = Cl/lj + C2/?2 + C3/I3, 

where 

Al = 1 (l _ 2 V3 + 2 2/3) > ^ = 1 (! + 2 l/3) ( 2 _ 2 l/3 ± 2 l/3 ^ 

are the roots of the characteristic polynomial z 3 -z 2 +z- 1 /3, and C ; are constants determined by the initial 
conditions. In the present case, as the proof of Lemma 2 shows, these initial values are given by 

11111 1 

Ci = -a? 1 = — , Co — —a-i 2 = — , c-i — —cka = . 

1 2 ' 4 1 3 ' 6 4 ' 12 

They can also easily be read off Table 1, where we have recorded the functions q^ for k = 0, . . .,4. The 

constants C, are thus the solution of the Vandermonde system 

q=Ci4+C 2 4+C34' i = 1,2,3, 



_,,„ " , 0<jc<1. (3.3) 
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Table 1. Values of the coefficients a^ , determining the finite-time distributions of the 
four-species Bak-Sneppen model. Values of /?;, i = 1,2,3, emphasised. Empty cells 
indicate zeros. 



that is 



C 3 +Cl/l2/l3-C2(/l2 + /l3) „ C3+CiAiA3-C 2 (Ai+A 3 ) c 3 + ci/tii 2 - c% (A\ + /t 2 ) 

c i = — m ttti — ' c 2 = tz rr-rr; — > <-3 



A l (A l -A 2 )(A l -A 3 ) 
It is then easy to compute q as 

00 3 



(A l -A 2 )A 2 (A 2 -A 3 ) 



(A 2 -Ai)A3(Ai - A3) 



xAj(2-xAj) 3 x(3-jc(3-jc)) 



,= 1 



.7=1 i=l 



(1-x^) 2 2(3-x(3-x(3-x))) 



2 ' 



Remark 1. The form of the function g with all Horner coefficients being equal to 3 is very interesting. We 
do not have an intuitive explanation for this peculiarity. 

We have thus proved the following result. 

Theorem 1. Let q be as in Proposition 4. The stationary distribution p — 1x4 of the Bak-Sneppen model 
with four species is given by 



,(d 4 *) = 



Z 



l [0,U 



{x)<Tx. 



(3.4) 



In particular, the one-dimensional margined of n, which is the asymptotic fitness distribution at a fixed site, 
is absolutely continuous with density 

(3.5) 



— P(/l,oo <*) = ^ +9(*)jl[ftl](jf). 

In Fig. 1 we have plotted the densities of the one-dimensional marginal distributions of f k for k = 0, . . . , 6, 
together with their limit as given by Eq. (3.5). The convergence asserted by Proposition 4 is clearly visible. 



4. Further steps towards a solution for general population sizes 



For general values of n, a similarly complete analysis of the asymptotic fitness distribution in the BS 
model as we presented for n — 4 in the previous section seems difficult to obtain. Here we collect some 
formulae that serve as the starting point if one wants to carry through the same approach as in Section 3. 
Similarly to Eq. (3.1), it holds for every v = !,...,« that 



Z a "j 



I I 

,/c+l,o-,i- c j y [ A {0<j: ]v[i[r( i ) <...<j: ]v[jO - ( „_3 ) ) 



hi J[°.l] 3 



l {f 2 <min(fi, ft ,*]„[)} 



{*Mf) 



(1,1, cr 



d 3 f 



(4.1) 
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Figure 1. Plot of the densities of the one-dimensional marginal distributions of / 4jt for 
k = 0,. ..,6 (dashed), together with their limit (solid) as given by Eq. (3.5) 



The integrals that appear on the right hand side of this equation fall in four categories, depending on how 
many of the three £'s are hit by the outer projection (•)]«[■ 

Case 1: // = 2. In this case, all three £'s are hit by the outer projection and the integral on the right side of 
Eq. (4.1) simplifies to 

f 3 1 {&<m«(ft^3^y[)} (*m)' 1 {0<i ]F[MI) <...<x ] „ l ^ ( „_ 3) <l} d3 f 



/[0,1] 3 
= (*]v[) - C ]v[,cr(l) 



3*]v[,<r(l) X ]v[,oiU + 1 



{ <*M,cK1) < -" ,Q: M,o-(b-3)< 1 } - 



Case 2: fx e {1,3}. In this case, either ft (if /j. = 3) or ft (if fi - 1) is not hit by the outer //-projection. For 
definiteness, we assume that fi — 3, the other sub-case being analogous. The sequential application of the 
outer projections (0]v[^ and (-)]3[ to x results in the vector y = (x v+ j, ...,x n ,X\,.. . ,Xy-2>ft), contributing a 



factor gf 3 to the integral. With the definitions y : 



cr(O) - : 



: 0, yaXn. 



(n-2) 



1 the integral in Eq. (4. 1) thus becomes 



■M:n-4 



f 

J[0,1] 3 



^{&<min(fl.6^]v[)} . 



y <r(o-- 1 (n-3)- 1) 1 ^o^tr- 1 (n-3)+l) 



(*) 

{0OV (l )C..«JV(„-3)«l) 



(4.2) 



where y l . n _ 4 = (x v+ ^,...,x„,xi,...,x v -2), iv.n-A = Oi, ■ ■ ■ ,'n-4), and the superscript (*) indicates that only 
the ordering of the components of y that pertain to x, in contrast to is restricted. Integration of this 
expression is standard and left to the reader. 

Case 3: // e {«,4}. In this case, both ft an d either £3 (if /j. — 4) or ft (if pi-ri) are hit by the outer 
projection. Again, for definiteness, we assume that fi — 4, and observe that the other sub-case can be 
treated in an analogous way. In this case 



( x M f ) 14| =y : = (Xy+4,...,X n ,XU---,Xy-2,{;i,fr), 



/]4[ 

resulting in a factor ft" -4 ^" -3 m the integrand; the integral in Eq. (4. 1) is clearly zero unless cr(l) = n - 3, 
so as not contradict the first indicator function that specifies ft to be less than ft and X\ V [. With the obvious 
definitions of Ji „_5 and i\-. n -5, the integral then becomes 



y »l:n-5 
•M:n-5 



Xa]3 lfe< ^ (ft ^ v[)}1 k 



(4.3) 



t(o-- 1 (n-4)-l) ^ sS - ) 'o-(o-- 1 (n-4)+ 1 ) / 

the evaluation of which we leave to the reader. 

Case 4: all other /j. In this last case, all three £'s remain unaffected by the //-projection and contribute to 
making the integral a little more complicated than before. Similar to the previous cases one obtains 

y : - (*Mf) M = [Xfi+2, ••• ,X„,X\, ... ,X v -2,£l, ft, ft, Xy+2, ... ,X M _ 2 ) , 



s 
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and, again, the integral vanishes if JVn ^ ft- Consequently, the integral in Eq. (4.1) becomes 



'•If 

IV[0,i] 



3 %2<min(£ 1 ,f 3 ,* M )} 1 ( 



{0<^ (2) ^...<3- [r( „_3 ) <l}' 



1 (v+ 1 )- i) 3 I (v+1)+1 ) ) 



(4.4) 



where = (x jU+2 ,...,x n ,jci,...,x v _2,Jc v +2,...,x jU _ 2 ), 3> A = & = 1,2,3, and i* = (ij)j*j k . According to 
whether or not %\ and £3 are separated by some xj in (y^n^j, there are a number of sub-cases to consider, 
which, however, do not present any difficulty. Closed-form expressions of the integrals 4.2, 4.3, and 4.4 
have been obtained, but are omitted for lack of space and because they do not add much insight at the 
current stage. 

By equating equal powers of x, Eq. (4.1), can thus, in principle, be used to derive an explicit recursion 
for the coefficients a n ^cr,i- Here, we only illustrate this potential by giving the recursion for n = 5, leaving 
the problem of working out the combinatorics of the general case for future research. An easy symmetry 
argument shows that Qf5 1 fe,(i,2),(i 1 ,i 2 ) is equal to a5X0-,\),(h,i\) so triat ^ i s enou gh to consider the identity 
permutation id e £2. 

Proposition 5. For n — 5, the coefficients akjj '■= tt5,*,id,i have the following properties: 

i) a k,0,j — Ofor all k > 1 and j > 0. 

ii) For k > 1 and j > 1, they satisfy the recursion 



V 3A-+1 "k.jj, 
t-*p=0 p+l 

a k,l,j~ 2 a k+U,j 



j-l 2p-;+l 



Jp=0 (p+V)(j-p) ak 'i- l -P>P> 
■j' 

\ptkj-ij - [l + ^=ij] ak,i-2,j + [5 + ^]a fc ,i- 3 j, 
Hi) For £ > 1 and j = 0, they satisfy ak,ifl -0 as well as the recursion 

(9 y.34+1 a k,\,p 

z 2j p= o p+i . 

O%,i-l,0 _ [l + (PI)/] a k,i-2fi + [3 + (P2)/] ff *,i-3,0 
(l X 2\ Vl-2 a kj-2-p,p , (I 2\ V i-3 "t.i-3-p.; 

"l i+ /;^p=o p+i + ^2 + /y^p=o p+i 

v i-2 a k,i-2-p.p 1 v ,'-3 a kj-3-p,p 



i = 1, 
z = 2, 
i>3. 



i = 2, 



(4.5) 



(4.6) 



i>3. 



Unfortunately, it does not seem possible to find an explicit solution to this recursion. In fact, the algo- 
rithm HYPER for linear recurrence relations with polynomial coefficients [15] indicates that the last line 
of Eq. (4.5) does not have a hypergeometric solution. Our results can, however, be used to obtain exact 
expression for the finite-time fitness distributions in the five-species BS model without having to evaluate 
any integrals. Calculations based on Eq. (4.5) and Eq. (4.6) give rise to the following conjecture. 

Conjecture. The coefficients a^u = a$j^i have the following properties: 
i) They form a stair case pattern in the sense that, for j > 0, 



ak,ifi 

a k,i,i 



i > 3k, 
i>3(k-l), 



0, i>3{k-A), A- 



7+1 



in particular, ak,\j = for j > 3k— 1. 
ii) They become constant at time 2 + jf+1, that is there exist numbers /?y such that arjyj = Pij for all 
k>i + 

In fact, we conjecture that, for each n, there exist k„ j|Tj ,- e N and p n ,cr,i e Q such that a n ,k,a-,i = Pn,a-,i for all 
k ^ k n (T i. 
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